Instability and Evolution of Nonlinearly Interacting Water Waves 
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We consider the modulational instability of nonlinearly interacting two-dimensional waves in deep 
water, which are described by a pair of two-dimensional coupled nonlinear Schrodinger equations. 
We derive a nonlinear dispersion relation. The latter is numerically analyzed to obtain the regions 
and the associated growth rates of the modulational instability. Furthermore, we follow the long 
term evolution of the latter by means of computer simulations of the governing nonlinear equations 
and demonstrate the formation of localized coherent wave envelopes. Our results should be useful 
for understanding the formation and nonlinear propagation characteristics of large amplitude freak 
waves in deep water. 
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Recently, there has been much interest 0, 0, S S 
in investigating the nonlinear formation of freak waves 
(also known as rogue waves, killer waves, or giant waves) 
in the oceans. Such waves are responsible for the loss 
of many ships and lives Q. Ocean waves usually have a 
spectrum of waves that depends on the variations in the 
wind speed and direction |S( • Freak waves are extraor- 
dinarily large amplitude localized water surface excita- 
tions, whose heights exceed many times the wavetrain 
height 9]. They may occur both in deep and shallow 
water, and may be created due to both statistical, linear 
(geometrical and spatio-temporal focusing) and nonlin- 
ear effects 6]. In nonlinear models, these steep objects 
are created on account of a delicate balance between the 
nonlinearity of the fluid and wave dispersion. The nonlin- 
ear dynamics of a single modulated large-amplitude wa- 
ter wavetrain can be modeled by a nonlinear Schrodinger 
equation 0,0,0,0,^^, which admits the Benjamin- 
Feir (modulational) instability [01 EE EDI EB] and the 
formation of rogue waves . Even though a wide spec- 
trum in the wave field may stabilize a wave train, the 
Benjamin-Feir instability is quite robust against a nar- 
row spectrum random field [6]; it occurs if the spectral 
width a of the wave is less than twice the average steep- 
ness (defined as k^a where ko is the central wavenumber 
and a is the rms value of the amplitude) [2(j. We stress 
that collapse of waves is inherently a two-dimensional 
problem, and that the single wave case has been inves- 
tigated extensively in Ref. E3- This indicates that the 
higher dimensionality of the problem at hand introduces 
a new flavor to the dynamics of interacting deep water 
waves, and brings us closer to observational prerequisites 
for such waves [2lj . 

rr 

Onorato et al. 3} recently developed a simple weakly 
nonlinear model for two nonlinearly interacting water 
waves in deep water with two different directions of prop- 
agation. They showed that the dynamics of these coupled 
waves is governed by two coupled nonlinear Schrodinger 



(CNLS) equations, and presented an investigation of the 
modulational instability for a one-dimensional two-wave 
system. A fourth-order nonlinear evolution equation for 
two Stoke wave trains in deep water was derived by Dhar 
and Das 22]. The two- wave system is similar to what 
has also been investigated in nonlinear optics [13, 124^ . 
in Bose-Einstein condensates [2^1 Hil, m signal trans- 
mission lines [22!, in nonlinear negative refraction in- 
dex meta-materials [^J |23, and in nonlinear plasmas 

In this Letter, we shall present for the first time a 
theory for the modulational instability of a pair of two- 
dimensional nonlinearly coupled water waves in deep wa- 
ter, as well as the formation and dynamics of localized 
freak wave packets. For this purpose, we shall use the 
CNLS equations of Onorato et al 3], which are valid 
for a system of obliquely propagating waves (crossing 
sea states). In Ref. 3] the x— axis has been denned 
as the middle (dichotome) between the two directions 
of propagation, viz. = (kA,x, ^A,y) = and 

k# = (kB, x , kB, y ) = (&, —I)] we shall assume that both k 
and / are larger than zero. The frequencies ujj of the 
two carrier waves (i.e. j = A, B) are related to the 
wavevectors kj by the deep water dispersion relation [^H 
ujj = y/g\kj\, where g denotes the gravitational accelera- 
tion. We note that the above hypothesis on the two wave 
propagation directions implies uoa = oob = J~9^^ where 



we have defined the wavenumber norm k = \/k 2 + I 2 . 
Moreover, we investigate the full dynamics of nonlin- 
early interacting deep water waves subjected to modu- 
lational/filamentation instabilities. It is found that ran- 
dom perturbations can grow to form inherently nonlin- 
ear water wave structures, the so called freak waves, 
through the nonlinear interaction between two coupled 
water waves. The latter should be of interest for explain- 
ing recent observations in water wave dynamics. 

Multiplying by i the system of two-dimensional CNLS 
equations (4) and (5) of Ref. 3], and correcting a misprint 
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(the coefficient in front of the d 2 B /dxdy term in Eq. (5) 
of Ref. 3), we have 
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where A and B are the amplitudes of the slowly 
varying wave envelopes such that the surface elevations 
(in meters) are given by tja = (l/2)A(r, t) exp(ikx + 
ily — iut)-\-c.c. and tjb = (l/2)5(r, t) exp(ikx — ily — 
iut) -\-c.c. where c.c. denotes complex conjugate. Here, 
C x = uok/2n 2 and C y = ujI/2k 2 are the group speed 
components along the x and y axis, respectively, and the 
group velocity dispersion coefficients are a = uo(2l 2 — 
/c 2 )/8k; 4 , p = u(2k 2 - / 2 )/8^ 4 and 7 = -3^//c/4^ 4 , and 
the nonlinear ity coefficients are given by £ = uon 2 /2 and 
C = u(k 5 -k 3 l 2 -3kl 4 -2k 4 K + 2k 2 l 2 K + 2l 4 K)/2K 2 (k-2K) 
(see Ref. 3]). It should be pointed out that Eq. (1) was 
first written down about forty years ago by Benney and 
Newell |3(| who gave formulae for all the coefficients ex- 
cept £ and C, which however can be found in Ref. |3(. 

It is easy to verify that the system of two-dimensional 
CNLS equations (QJ has the space independent har- 
monic solutions A eq = A ex.p[—i(^\Ao\ 2 + 2(\B \ 2 )t] 
and B eq = B exp[-i(£\B \ 2 + 2(\A \ 2 )t}. Assum- 
ing a small (linear) harmonic perturbation around the 
above mentioned equilibrium states with the wavevec- 
tor K = (K, L) and the frequency i.e. substituting 
with A — > (A + eAi)exp[-z(£|A | 2 + 2(\B \ 2 )t] and 
B (B + eB 1 )exp[-i(^\B \ 2 + 2(\A \ 2 )t] into (1), lin- 
earizing in the small real parameter e <C 1, then separat- 
ing the real and imaginary parts, combining the resultant 
equations, and Fourier transforming, we obtain the non- 
linear dispersion relation 

[(n-c x K-c y L) 2 -n 2 ] [(n-c x K+c y L) 2 -n 2 ] =n 4 c , 

(2) 

where Qj = (aK 2 + PL 2 - ~/KL)(aK 2 + PL 2 + ^KL + 
2£|A | 2 ), n| = (aK 2 + /3L 2 + ^KL)(aK 2 + /3L 2 - 
7 KL + 2£ 1 5 1 2 ) , and fi* = 1 6 C 2 | A \ 2 \ B \ 2 (aK 2 + PL 2 - 
r yKL)(aK 2 + PL 2 + jKL). For one-dimensional wave 
propagation, i.e. for L = 0, Eq. is identical to Eq. 
(11) of Ref. p|. Onorato et al. [3 presented numerical re- 
sults for the modulational instability regimes and growth 
rates for L = and Ao — Bo from their Eq. (11). We 
note that the dispersion relation (2) is nonlinear in the 
wave amplitudes (but linear in the small expansion pa- 
rameter). 

Taking into account that a measure of the wave non- 
linearities are given by the "steepnesses" kA and kB, 



we will use the natural normalizations Ao = A' /k and 
Bo = B'o/n for the wave amplitudes. For the wave num- 
bers and frequencies we use K' = K/k, V = L/k, 
k' = k/n, V = 1/k, and fl f = Q,/uj. The coefficients are 
normalized as C' x = C x k/cu = fc , /2, C' y = CyKjuj — Z'/2, 
a' = an 2 /u = (2l' 2 - k ,2 )/8, p' = Pk 2 /u = (2k' 2 - / /2 )/8, 
i = >jk 2 /uj = Sl'k'/A, f = ijujk 2 = 1/2, and 
C = C/ook 2 = [(k') 5 - (k') 3 (l') 2 - 3fc'(Z') 4 - 2(^) 4 + 
2(k') 2 (V) 2 + 2(V) A )/2(k' - 2), and all quantities in Eqs. 
(2) and (3) can be replaced with their primed counter- 
parts. Here we have k! = cos 9 and V = sin#, where is 
the angle between the wave directions and the dichotome. 
Except for numerical factors, Eqs. (2) and (3) then only 
contain the angle the known wavenumbers K' and L', 
the wave amplitudes A' and B^ and the unknown fre- 
quency Q! . 

In the following, we numerically solve our nonlinear 
dispersion relation and present the growth rate Y (the 
imaginary part of fi) in Figs. 1-3, where we have studied 
the impact of different angles 6 on the growth rates for 
interacting waves. For all cases we used the waves am- 
plitudes 0.1/ft for the two waves. In the left-hand panels 
of Figs. 1 and 2 we show the single wave cases, which 
exhibit the standard Benjamin- Feir instability [cf. Eq. 
(3.6) and Fig. 1 of Ref. Q], tilted by the angle in the 
(K,L) plane. The right-hand panels show the cases of 
interacting waves. We see from Fig. 1 that a relatively 
small = 7r/8 gives rise to a new instability with a max- 
imum growth rate that is more than twice as large as 
the ones for the single wave cases, in the direction of the 
dichotome. For a larger angle # = 7r/4^0.79, displayed 
in Fig. 2, we see that the two waves do not interact to 
enhance the linear growth rate significantly. We note (as 
pointed out in Ref. that the coefficient a changes 
sign when 6 = arctan(//fc) = arctan(l/\/2) ~ 0.615 rad 
~ 35.3° so that we have a focusing (defocusing) insta- 
bility along the x axis for < 0.615 (> 0.615). This 
clearly stabilizes the waves so that they only exhibit the 
standard Benjamin- Feir instability but do not interact 
to enhance the instability. On the other hand, for two 
counter-propagating waves with = 7r/2, we see in Fig. 
3 that there is again an instability along the dichotome 
(in agreement with Ref. 0) and also obliquely in narrow 
bands almost perpendicular to the two waves. 

In order to study the dynamics of nonlinearly inter- 
acting water wave packets, we solve the coupled system 
of equations Q. The results are displayed in Fig. 4. In 
the numerical simulation, we have used the normaliza- 
tion A! — A/k, B f = B/k, t f = out, x' = kx, and y' = ny 
(the other scaled parameters are as above), while the re- 
sults are shown in dimensional units in Fig. 4. We have 
used the same parameters as in the right-hand panel of 
Fig. 1, where the two interacting waves initially have 
the amplitude A = B = 0.1/ft, plus low-amplitude noise 
(random numbers) of the order 10 -3 /ft to give a seed 
for any instability. The wavefronts in Fig. 4 are initially 
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FIG. 1: The normalized growth rate V /u) plotted as a function 
of K/k and L/k. Here we have used = tt/8 and the wave 
amplitudes Ao = Bo = 0.1/ k. The left upper and lower panels 
show the cases with a single wave Ao and Bo, respectively, 
while the right panel shows the case of interacting waves. 



directed primarily in the x direction, reflecting the max- 
imum growth rate in the K direction in the right-hand 
panel of Fig. 1. The linear growth saturates in the forma- 
tion of wavepackets, localized in the x direction, that are 
strongly correlated between wave A and B; see the pan- 
els at time t = 900/a;. At later times the waves A and B 
decouple and there are scattered large-amplitude waves 
that are localized both in the x and y directions. In order 
to compare with the case of a single wave we have also 
run a simulation where B was set to zero, so that we have 
the standard Benjamin- Feir instability shown in the up- 
per left panel of Fig 1. In this case (not shown here), we 
could not see the formation of well-defined wave-packets, 
but the instability gave rise to a wide spectrum of waves 
in different directions, in agreement with the linear anal- 
ysis in the left-hand panels of Fig. 1. The new instability 
due to the coupling of the two waves, shown in the right- 
hand panel of Fig. 1, has a well-defined maximum in 
the x direction, which is important for the localization 
of wave energy into localized wavepackets seen in Fig. 4. 
Taking some typical data from ocean waves Q where a 
typical wave frequency is 0.09 Hz, we have uo = 0.56 s -1 , 
and k = uj 2 / g ~ 0.033 m _1 . The waves A and B in Fig. 
4, then have the initial amplitudes \A\ = \B\ = 0.1/ k ^ 3 
meters. In Fig. 4, we see at t = 1200/a; (~ 670 seconds) 
that wave A has some localized wave packets with a max- 
imum amplitude of « 0.35/ft w 10 meters. 

To summarize, we have presented a theoretical study 
of the modulational instabilities of a pair of nonlinearly 
interacting two-dimensional waves in deep water, and 
have shown that the full dynamics of these interacting 
waves gives rise to localized large-amplitude wavepack- 
ets. Starting from the CNLS equations of Onorato et al. 
0, we have derived a nonlinear dispersion equation. The 
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FIG. 2: The normalized growth rate T/uj plotted as a function 
of K/k and L/k. Here we have used 6 — tt/4 and the wave 
amplitudes Ao = Bo = 0.1/ k. The left upper and lower panels 
show the cases with a single wave Ao and Bo, respectively, 
while the right panel shows the case of interacting waves. 
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FIG. 3: The normalized growth rate F/oo plotted as a function 
of K/k and L/k. Here we have used = tv/2 and the wave 
amplitudes Ao = Bo = 0.1/ k. The left upper and lower panels 
show the cases with a single wave Ao and Bo, respectively, 
while the right panel shows the case of interacting waves. 



latter has been numerically analyzed to show the depen- 
dence of obliqueness of the interacting waves and of the 
modulations on new classes of modulational instabilities 
that we found in multi-dimensional situations. Further- 
more, the numerical analysis of the full dynamical system 
reveals that even waves that are separately modulation- 
ally stable can, when nonlinear interactions are taken 
into account, give rise to novel behavior such as the for- 
mation of large-amplitude coherent wave packets with 
amplitudes more than three times the ones of the initial 
waves. This behavior is very different from that of a sin- 
gle wave which experiences the standard Benjamin- Feir 
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FIG. 4: The interaction between two waves, initially with equal amplitudes \A\ — \B\ — 0.1ft -1 and a propagation angle of 
= 7r/8 relative to the dichotome. Added to the initially homogeneous wave envelopes is a low-amplitude noise of order 10 -3 / k 
to give a seed to the modulational instability. 



instability which dissolves into a wide spectrum of waves. 
These results will have relevance to the nonlinear insta- 
bility of colliding water waves, which may interact non- 
linearly in a constructive way to produce large-amplitude 
freak waves in the oceans. 
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